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Abstract. In this talk I have presented the data analysis results of extracting properties of 
halo WIMPs: the mass and the (ratios between the) spin-independent and spin-dependent 
couplings/cross sections on nucleons by the AMIDAS website by taking into account possible 
unrejected background events in the analyzed data sets. Although non-standard astronomical 
setup has been used to generate pseudodata sets for our analyses, it has been found that, 
without prior information/assumption about the local density and velocity distribution of halo 
Dark Matter, these WIMP properties have been reconstructed with ~ 2% to < 30% deviations 
from the input values. 



1. Introduction 

In order to extract properties of halo WIMPs (Weakly Interacting Massive Particles) by using 
data from direct Dark Matter detection experiments as model-independently as possible, we 
have developed a series of data analysis method for reconstructing the one-dimensional WIMP 
velocity distribution [1] as well as determining the WIMP mass [2], the spin-independent (SI) 
WIMP coupling on nucleons [3] and the ratios between different WIMP couplings / cross sections 
[4]. Moreover, in collaboration with the DAMNED (DArk Matter Network Exclusion Diagram) 
Dark Matter online tool [5], part of the ILIAS Project [6], the "AMIDAS" (A Model-Independent 
Data Analysis System) website for online simulation/data analysis has also been established 
[7,8,9]. 

In this article, in order to demonstrate the usefulness and powerfulness as well as the model- 
independence of the AMIDAS package for direct Dark Matter detection experiments, I will analyze 
blindly some pseudodata sets generated theoretically for different detector materials and present 
the reconstructed WIMP properties. This means that I will simply upload these data sets onto 
the AMIDAS website and follow the instructions to reconstruct different WIMP properties without 
using any information about the input setup used for generating the analyzed pseudodata. For 
cases in which some information about WIMPs (e.g., the mass m x ) and/or Galactic halo (e.g., 
the local Dark Matter density po) is required, I will naively use the commonly used/favorite 
values for the data analyses. Moreover, due to the time limit on preparing the presentation for 
the CYNGUS 2011 workshop, in Ref. [10] we considered only the case of null background events. 
In this article, we have taken into account a small fraction of possible unrejected background 
events in the analyzed data sets [11, 12, 13, 14]. 
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Figure 1. The reconstructed WIMP mass m x by means of the method introduced in Ref. [2] 
with a target combination of 28 Si + 76 Ge nuclei. Two forms of the elastic nuclear form factor 
given in Eqs. (1) and (2) have been used in the left and right frames, respectively. 



After that I show the blindly reconstructed properties of halo WIMPs in Sec. 2, in Sec. 3 I 
will reveal the input setup and the background spectrum used for generating the analyzed data 
and compare the reconstructed results to them. Finally I conclude in Sec. 4. 



2. Reconstructed WIMP properties 

In this section, I present the reconstructed WIMP properties analyzed by the AMIDAS website. 
While in each uploaded file there are exactly 50 data sets, in each data set there are on average 
50 recorded events (i.e., 50 measured recoil energies); the exact number of total events is Poisson 
distributed. For simplicity, the experimental minimal and maximal cut-off energies have been 
set as and 100 keV for all data sets. 

In order to check the effect of using a "wrong" elastic nuclear form factor, two forms have 
been considered for the SI WIMP-nucleus cross section in our analyses. One is the simple 
exponential form: 

^ex(Q) = e" Q/Q ° • (1) 
Here Q is the recoil energy transferred from the incident WIMP to the target nucleus, Qq is the 



fm 



nuclear coherence energy given by Qo = 1.5/ton-Ro, where Rq = 0.3 + 0.91 (m-^/GeY) 1 ^ 3 
is the radius of the nucleus and tbn is the mass of the target nucleus. Meanwhile, we used also 
a more realistic analytic form for the elastic nuclear form factor: 
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Here j\(x) is a spherical Bessel function, q = y/2m^Q is the transferred 3-momentum, for the 
effective nuclear radius we use 



Rt = ^R\ - 5s 2 (3) 

with Ra — 1-2 A 1 / 3 fm and a nuclear skin thickness s ~ 1 fm. For the SD WIMP-nucleus cross 
section, we only used the "thin-shell" nuclear form factor: 

pi (q) = { 3o(lRl) > f ° r QRi < 2 -55 or qR Y > 4.5 , 

Tsl ' } const. ~ 0.047 , for 2.55 < qRi < 4.5 . U 



2.1. WIMP mass m x 

As one of the most important properties of halo WIMPs as well as the basic information for 
reconstructing other quantities in our model-independent analysis methods, I consider at first 
the determination of the WIMP mass m x by means of the method introduced in Ref. [2]. 

In Figs. 1 I show the reconstructed WIMP masses and the lower and upper bounds of their 
la statistical uncertainties. The usual target combination of 28 Si + 76 Ge nuclei has been used 




Figure 2. The reconstructed squared SI WIMP-nucleon coupling |/ p | 2 by means of the method 
introduced in Ref. [3] with a 76 Ge target. The commonly used value of the local Dark Matter 
density po = 0.3 GeV/cm 3 and a larger value of p$ = 0.4 GeV/cm 3 as well as the elastic nuclear 
form factors given in Eqs. (1) and (2) have been used for estimating \ f p \ 2 in the left and right 
frames, respectively. 



for this reconstruction, whereas two analytic forms of the elastic nuclear form factor given in 
Eqs. (1) and (2) have been used for determining m x in the left and right frames, respectively. 
While m x , n with n = —1, 1, 2 and m Xi<7 have been estimated by Eqs. (34) and (40) of Ref. [2], 
respectively, ?n x ,combined has been estimated by the x 2 ~fitting defined in Eq. (51) of Ref. [2], 
which combines the estimators for m Xjn and m Xt(T with each other. The reconstructed WIMP 
mass m XjCom bined as w ell as m Xi „ and m XtCT shown here have been corrected by the iterative 
<3max~matching procedure described in Ref. [2]. 

It can be found here that, although all single estimators (m x ^ n with n = —1, 1, 2 and rn Xt0 ) 
give generally a (relatively lighter) WIMP mass of ~ 60 GeV or even lighter and a lex upper 
bound of ~ 140 GeV, the mean values of the combined (in principle, more reliable) results 
(second columns) of the reconstructed WIMP mass give m x ~ 115 GeV with a rough la upper 
(lower) bound of ~ 190 (80) GeV, or, equivalently, 

m x ~ 115+H GeV . (5) 

Moreover, the combined results with two different form factors show not only a large overlap 
between ~ 85 GeV and ~ 180 GeV, but also a good coincidence: comparing to the ~+35 GeV 
la statistical uncertainty and the 

^^30 GeV overlap, the difference between two median values 
is < 30 GeV! This indicates that, for the first approximation of giving/constraining the most 
plausible range of the WIMP mass with pretty few total events, the uncertainty on the nuclear 
form factor could be safely neglected. 

2.2. Spin-independent WIMP-nucleon coupling |/ p | 2 

Following the WIMP mass determination, I consider now the reconstruction of the SI WIMP 
coupling on nucleons |/ p | 2 with a 76 Ge target [3] 1 . 

In Figs. 2 I show the reconstructed squared SI WIMP-nucleon couplings and the lower and 
upper bounds of their la statistical uncertainties estimated by Eqs. (17) and (18) of Ref. [3] 
with an assumed (100±10 GeV, labeled with the subscript "input") and the reconstructed (from 
Sec. 2.1, labeled with "recon") WIMP masses. The commonly used value of the local Dark 
Matter density p$ = 0-3 GeV/cm 3 and a larger value of po = 0.4 GeV/cm 3 [15, 16, 17] as well 
as the elastic nuclear form factors given in Eqs. (1) and (2) have been used for estimating |/ p | 2 
in the left and right frames, respectively. 

Among these results, the mean value and the overlap of two most plausible results (estimated 

1 Remind that the theoretical prediction by most supersymmetric models that the SI scaler WIMP couplings on 
protons and on neutrons are (approximately) equal: f p ~ / n has been adopted in the AMIDAS package. 
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Figure 3. The reconstructed ratio between two SD WIMP-nucleon couplings, a n /a p , by means 
of two methods introduced in Ref. [4]. As usual, the elastic nuclear form factors given in Eqs. (1) 
and (2) have been used for determining a n /a p in the left and right frames, respectively. 

with the reconstructed WIMP mass) give a rough la range of 

|/ p | 2 ^ 1.05^;| xl(T 17 GeV- 4 , (6) 

or, equivalently, 

|/ p | ~ 3.24±g;|^ x 1(T 9 GeV~ 2 . (7) 

Since the reconstructed WIMP mass given in Sec. 2.1 is m x ~ 115 GeV, one can simply use the 
proton mass m p to approximate the WIMP-proton reduced mass m riP and give a reconstructed 
SI WIMP-nucleon cross section as [10] 

4p = (£) <M 2 « (£) < I/pI 2 - 5.03±i;f| x 10- 9 pb. (8) 

2.3. Ratio of two spin-dependent WIMP-nucleon couplings a Q /a p 

In Figs. 3 I show the reconstructed a n /a p ratios and the lower and upper bounds of their la 
statistical uncertainties estimated by Eqs. (2.7) and (2.12) of Ref. [4] with n = 1 as well as by 
Eqs. (3.16) and (3.20) of Ref. [4] at the shifted energy points [1, 4]. A combination of 19 F + 127 I 
targets has been used for the reconstruction of a n /a p under the assumption that the SD WIMP- 
nucleus interaction dominates over the SI one (labeled with the superscript "SD"), whereas a 
third target of 28 Si has been combined with 19 F and 127 I for the case of the general combination 
of both SI and SD WIMP interactions (labeled with the superscript "SI + SD"). 

It can be found that, firstly, the "+ (plus)" solutions of the a n /a p ratios given here are 
obviously too large to be the reasonable choice for a n /a p and the "— (minus)" solutions should 
be the correct ones 2 . Secondly, although the reconstructed results under the assumption of the 
SD dominant WIMP interaction (third columns) is in general larger than the (in principle more 
plausible) results obtained without such a prior assumption (last columns), one could still use 
the mean value and the overlap of these two results to give a rough la range of 

? - 0.65«1 . (9) 

Up 

2.4- Ratios of the SD and SI WIMP-nucleon couplings a x(pn)/°xp 

In Figs. 4 I show the reconstructed a ^ pn )/ a xp ra tios and the lower and upper bounds of their 
la statistical uncertainties estimated by Eqs. (3.9), (3.10) and (3.21) of Ref. [4] (with a n /a p 
estimated by Eq. (3.16) of Ref. [4]) as well as by Eqs. (3.25) and (3.29) of Ref. [4] at the 
shifted energy points. A combination of data sets of 19 F, 127 I and 28 Si targets (labeled with 

2 Remind that, as discussed in Ref. [4], the correct choice from the "+" and "— " solutions could be decided 
directly by the values of the group spins of protons and neutrons of the used target nuclei, (S( p ^). 
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Reconstructed ratio between the SD and SI WEVIP-neutron cross sections 
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Reconstructed ratio between the SD and SI WEVIP-neutron cross sections 
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Figure 4. The reconstructed ratios between the SD and SI WIMP-nucleon couplings, 
°x(p n) / °xp ' means of two methods introduced in Ref . [4] . As usual, the elastic nuclear 
form factors given in Eqs. (1) and (2) have been used for determining cr^P n )/c^p in the left and 
right frames, respectively. 



the superscript "XYZ") and that of data sets of 23 Na or 131 Xe with the (common) one of 76 Ge 
(labeled with the superscript "XY") have been used and the mean value and the overlap of these 
two results give a rough la range of 

% - 8.94±^ 3 x 10 5 , % ~ 3.16^-36 x 10 5 . (10) 

Then, firstly, from these results one can further obtain that [10] 3 

^ 0.59±g;?I . (11) 



Secondly, combining the results in Eq. (10) with given in Eq. (8), one can obtain that [10] 

4° ~ 4.50l};^ x 1(T 3 pb , 4° ~ 1.59±J;gl x 10~ 3 pb . (12) 

These results give in turn that [10] 

|a p | = 0.112±Hi, K\ = 0.067™. (13) 

On the other hand, one can also use the reconstructed a n /a p ratio given in Eq. (9) and one of 
the two results given in Eq. (12) to obtain that [10] 

a™ ~ 3.76± 3 J? x 1(T 3 pb , ~ 1.90±^ x 1(T 3 pb . (14) 

These results can also give that 

lopi- o.io3i8:8», Ki- o.o73t8:8^8 • (is) 

It can be found that, not surprisingly, the statistical uncertainties on the reconstructed cr^P, n ^ 
given in Eq. (14) are ~ 2 or 3 times larger than those given in Eq. (12): Since CT x(pn)/°xp 

3 Remind that the results given in the second and third columns of the tables in Figs. 4 are reconstructed with 
the On/dp ratio given in the last columns of the tables in Figs. 3. 



reconstructed with the F + I + Si combination involve already the reconstructed a Q /a p ratio 
given in Eq. (9), the uncertainties on ^(pn) gi ven m Eq. (14) are thus overestimated. Secondly, 
the reconstructed a n /a p ratio given in Eqs. (9) and (11) and the reconstructed ^(p 11 ) & ven * n 
Eqs. (12) and (14) seem to match to each other pretty well. 

The analyses given here show that, firstly, once one can estimate the SI WIMP- 
nucleon coupling/cross section, |/ p | or a^ p , and (one of) the ratios between the SD and SI 
WIMP-nucleon cross sections, and/or the ratio between two SD WIMP-nucleon couplings, 
the other couplings/cross sections could in principle be estimated. Secondly, the WIMP 
couplings/cross sections estimated in different way would be self-cross-checks to each other 
and the (in) compatibility between the reconstructed results would also help us to check the 
usefulness of the analyzed data sets offered from different experiments with different detector 
materials [14]. 



3. Input setup for generating pseudodata 

In Table 1 I give finally the input setup for generating the pseudodata sets used in the analyses 
demonstrated in the previous section. For comparison, the reconstructed results shown in the 
previous section are also summarized here. 

For generating WIMP signals, we used the commonly used form for the nuclear form factor 
given in Eq. (2) with another often used analytic form for R±: 
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with Ra ~ (1.23 A 1 / 3 - 0.6) fm, r ~ 0.52 fm, s ~ 0.9 fm. Moreover, the shifted Maxwellian 
velocity distribution: 
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with the Sun's Galactic orbital velocity vq = 230 km/s has been used; v e is the time-dependent 
Earth's velocity in the Galactic frame: 



v e {t) = v 



1.05 + 0.07 cos 



i yr 



(18) 



where the date on which the Earth's velocity relative to the WIMP halo is maximal has been set 
as t p = 140 d. In addition, different from our setup used in Ref. [10], the experimental running 
date has been set as t cxp t = 100 d and thus v e (t exp t) = 254 km/s, much larger than the usually 
used values: 200 km/s < v e < 240 km/s. Although these values for the astronomical setup 
are non-standard, we would like to stress that, as shown in the previous section and Table I, 
such a non-standard halo (model) would not affect the reconstructed results, since for using 
the AMIDAS package and website to analyze (real) data sets, one needs only the form factors for 
SI and/or SD WIMP-nucleaus cross sections, prior knowledge/assumptions about the WIMP 
velocity distribution fi(v) and local density po (except the estimation of the SI WIMP-nucleon 
coupling |/ p | 2 ) are not required. 

On the other hand, for generating background events, the target-dependent exponential form 
for the residue background spectrum introduced in Ref. [11]: 



dR\ 

dQ J bg,ex 



exp 



Q/keV \ 
^40.6 J 



(19) 
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Table 1. The input setup for generating the pseudodata sets used in the analyses demonstrated 
in this article. The theoretically estimated values and the reconstructed results are also given, 
f: estimated for 130 GeV m x ; X- different from Ref. [10]. 



has been used. Here A is the atomic mass number of the target nucleus; the power index of 
A, 0.6, is an empirical constant, which has been chosen so that the exponential background 
spectrum is somehow similar to, but still different from the expected recoil spectrum of the 
target nuclei (see Figs. 5) 4 . Meanwhile, considering possible radioactivity with a characteristic 
energy, we combined the exponential background spectrum in Eq. (19) with a Gaussian excess: 



dQJ bg.Gau 



2vr(a Q , bg /keV) 



exp 



(Q Qbg,peak) 



(20) 



4 Note that, among different possible choices the atomic mass number A has been chosen as the simplest, unique 
characteristic parameter in the general analytic form (19) for defining the residue background spectrum for 
different target nuclei. However, it does not mean that the (superposition of the real) background spectra would 
depend simply /primarily on A or on the mass of the target nucleus, ttin. 
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Figure 5. Measured energy spectra (solid red histograms) for a 19 F (left) and a Ge (right) 
targets. The dotted blue curves are the elastic WIMP-nucleus scattering spectra, whereas 
the dashed green curves are the exponential background spectra combined with the Gaussian 
excesses (rcau = 1> Qbg,peak = 50 keV and o"Q,b g = 1 keV) and normalized to fit to the chosen 
background ratio = 12%. See the text for further details. 



where Qbg,peak and 0Q,bg are the characteristic energy and energy dispersion of this background 
excess, respectively; rcau is the ratio of this Gaussian excess to the exponential spectrum. 

In Figs. 5 I show the measured energy spectra (solid red histograms) for a 19 F (left) and a 
76 Ge (right) targets. The dotted blue curves are the elastic WIMP-nucleus scattering spectra, 
whereas the dashed green curves are the exponential background spectra given in Eq. (19) 
combined with the Gaussian excesses given in Eq. (20) (rQau = 1> Qbg,peak = 50 keV and 
a Q,bg = 1 keV for both targets), which have been normalized so that the ratios of the areas 
under the background spectra to those under the (dotted blue) WIMP scattering spectra are 
equal to the background-signal ratio (rb g = 12%) in the whole data sets. The experimental 
threshold energies have been assumed to be negligible and the maximal cut-off energies are set 
as 100 keV. 500 total events on average in one experiment have been simulated. 

It can be seen here that 5 , firstly, due to the contribution of the unrejected background events 
in low energy ranges (Q < 20 keV for 19 F and Q < 40 keV for 76 Ge), the counting rates at 
the first energy bins (one of the two important quantities required in our model-independent 
data analyses) have been (strongly) overestimated. Secondly, the Gaussian background excesses 
around Q = 50 keV cause clearly overestimates of the event rates, which would not only 

contribute (significantly) to the estimates of /„ = X^aii events I F 2 (Qi) [1], but also cause 

a larger statistical fluctuation [2] . 

Nevertheless, our results obtained by analyzing pseudodata sets of O(50) total events showed 
that, firstly, the WIMP mass given in Eq. (5) can match the input value very well: the deviations 
between the input and the reconstructed values with different assumed nuclear form factors are 
< 20% ( < 30 GeV). As discussed earlier, this indicates that, for the first approximation of 
giving/constraining the most plausible range of the WIMP mass with pretty few total events, 
the uncertainty on the nuclear form factor could be safely neglected. 

Secondly, all WIMP-nucleon couplings/cross sections as well as the ratios between them have 
also been reconstructed with ~ 2% to < 30% deviations from the input /theoretically estimated 
values. Although the SI WIMP coupling |/ p | estimated with the larger (input) local Dark Matter 
density (right frame of Figs. 2) could be underestimated [3], one can at least give an upper bound 



5 More detailed illustrations and discussions about the effects of residue background events on the measured 
energy spectrum can be found in Refs. [11, 14]. 



on |/ p |. Moreover, by combining different methods for estimating different (ratios between the) 
WIMP couplings/cross sections, one could in principle observe/confirm the (in) compatibility 
between these results and probably correct the reconstructed values. 

4. Summary 

In this article I demonstrated the data analysis results of extracted WIMP properties by 
using theoretically generated pseudodata for different target nuclei, taking into account some 
unrejected background events mixed in the analyzed data sets. As an extension as well as 
the complementarity of our earlier theoretical works, I combined reconstructed results of the 
(ratios between different) WIMP couplings/cross sections on nucleons to estimate each individual 
coupling/cross section. Hopefully, the AMIDAS package and website as well as this demonstration 
can help our experimental colleagues to analyze their real direct detection data in the near future 
and to determine (at least rough ranges of) properties of halo Dark Matter particles. 
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